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This paper deals with the dynamical system that generalizes the MIC-Kepler system. It is shown that 
the Schrodinger equation for this generalized MIC-Kepler system can be separated in prolate spheroidal 
coordinates. The coefficients of the interbasis expansions between three bases (spherical, parabolic and 
spheroidal) are studied in detail. It is found that the coefficients for this expansion of the parabolic basis 
in terms of the spherical basis, and vice-versa, can be expresses through the Clebsch-Gordan coefficients 
£f~^ ' for the group SU(2) analytically continued to real values of their arguments. The coefficients for the 

£\j i expansions of the prolate spheroidal basis in terms of the spherical and parabolic bases are proved to 

satisfy three-term recursion relations. 



1 Introduction 

The generalized MIC-Kepler system is described by the equation £Q 

; 2 1 , ci c 2 



\ (-iV - sAf $ - 



where c\ and c 2 nonnegative constants, and 



i> = Ei>, (1.1) 



ctf 1 

3. A = - dv,~x,0), and rotA=— . (1.2) 

r(r — z) r° 
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(We use the system of units for which % = m = e = c— 1.) The monopole number s satisfies the Dirac's rule 
of charge quantization s — 0,±1/2,±1, . . .. Each value of s describes its particular generalized MIC-Kepler 
system. The Schrodinger equation l|l.l|l for Cj = Q(i = 1,2) and s ^ reduces to the Schrodinger equation 
of the MIC-Kepler system The MIC-Kepler system could be constructed by the reduction of the 
four-dimensional isotropic oscillator by the use of the so-called Kustaanheimo-Stiefel transformation both 
on classical and quantum mechanical levels [3] . In the similar way, reducing the two- and eight- dimensional 
isotropic oscillator, one can obtain the two- [5| and five-dimensional [Hj analogs of MIC-Kepler system. An 
infinitely thin solenoid providing the system by the spin 1/2, plays the role of monopole in two-dimensional 
case, whereas in the five-dimensional case this role is performed by the SU (2) Yang monopole 7 , endowing 
the system by the isospin. All the above-mentioned systems have Coulomb symmetries and are solved in 
spherical and parabolic coordinates both in discrete and continuous parts of energy spectra [HIE]- There are 
generalizations of MIC-Kepler systems on three-dimensional sphere an d hyperboloid as well. The 
MIC-Kepler system has been worked out from different points of view in Refs. |12l 1131 ITU 1151 llfi| . 

At s = 0, Eq. i|l.l|) is reduced to the Schrodinger equation for generalized Kepler-Coulomb system |17j . 
In the case when s = and ei = c 2 ^ 0, the equation (|l.lf> reduces to the Hartmann system that has been 
used for describing axially symmetric systems like ring-shaped molecules 

The system described by the Schrodinger equation 1)1. Q is one of the Smorodinsky-Winternitz type 
potentials |19j . The Smorodinsky-Winternitz type potentials where revived and investigated in the 1990 by 
Evans [20]. 

In Ref. PP it is shown that the variables in Schrodinger equation i|l.l|) are separated in spherical and 
parabolic coordinates. In this article is shown that the variables in Eq. Ijl-ljl can be separated in the prolate 



1 



spheroidal coordinates also. The system of spheroidal coordinates a natural system for investigating many 
problems in mathematical physics (see Ref. |21| and references cited therein). In quantum mechanics, the 
spheroidal coordinates play an important role because they are appropriate in describing the behavior of 
a charged particle in the field of two Coulomb centers. The distance R between the centers is a dimen- 
sional parameter characterizing the spheroidal coordinates, these coordinates are changed into spherical and 
parabolic coordinates as R — > and R — ► oo respectively, if the positions of one Coulomb center and the 
charged particle are fixed when taking the limits. In this sense, the spheroidal coordinates are more general 
then the spherical and parabolic coordinates. 



2 Spherical and Parabolic Bases 

For completeness, we here present the solutions of the Schrodinger equation found in pQ. Eq. in 
the spherical coordinates becomes 



d 2 



4 cos 2 % \dip 



4ci 



4 sin 



2 



— + 2is 



4c 2 



2[E+-}}i/j = Q, 



(2.1) 



where 



A,. 



1 d ( d 



1 d 



sin 6 



r 2 dr V dr J ' sin dO \ 86 
The solution of Eq. (|2.1I has the form 

Vfil (r, 0, <p; 6 U 6 2 ) = R% (r; 6 U 6 2 ) Z$ (0, tp; 6 U 6 2 ) . 

The functions (9, cp; 61,62) and R^) (r; 61,62) are givenly the formulae 



(2.2) 



(2.3) 



As) 1 



),cp;6i,6 2 ) = N jm (6 1 ,6 2 ) 



cos - 
2 



sin - 
2 



(2.4) 



R ( n](r) = C nj (6 1 ,6 2 )(2ery+^e- Er F (-n + j + 1; 2j + S 1 + 6 2 + 2; 2er) , (2.5) 

where Pn^\x) are the Jacobi polynomials, F(a;c;x) is the confluent hypergeometric function, Nj m {5\,5 2 ) 
and C n j(6i,62) are normalization constants 



AT ( A A x (2j + 61 + 62 + - m+)\T(j + m+ + Si + 62 + 1) . . 

1/ 47iT(j - m_ + 6i + l)T(j + m_ + 6 2 + 1) 



2e 2 T(n + j +8 1 + 5 2 + l) 

r(2j + 6i + 6 2 + 2)V 



We assume that 

As) fn r x„\ y( s ) 



BinOZ)^, {e,<p;8 u 82)Z)Z (0, tp; 8 U 6 2 ) dO d V = 8 3f 5 mm ,, (2.. 



r 2 (r; ^ , 6 2 ) R[ s ] (r; ft, fc) dr = <f„ 



(2.9) 
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and denote by the e the following expression 



The energy spectrum has the form 



e = V^2E = 



E = EP = -- 



and the quantum numbers to and j run through the values: m = —j, — j + 1, . . . , j — 1, j and 



(2.10) 



(2.11) 



J 



\m + s\ + \m — s \m + s\ + \m — s 



2 ' 2 

We make the following notation also m± = (|m + s\ ± \m — s|)/2 and 



1. 



To! = |m - s\ + Si = \J (m - s) 2 + 4ci, m 2 = |m + s| + S 2 = \J (to + s) 2 + 4c 2 
The wave functions (|2.3() are the eigenfunctions of commuting operators M and J z , moreover and 



M$L(r t 0,<p-,8 1 ,5 2 )= lj + 



<5i + <5 2 



where 



M = J z 



2ci 



3 + ^4^ + l)$UrA<p-,Si,S a ), 



2c 2 



1 + cos 9 1 — cos ( 
Here J 2 is the square of the angular momentum 

i 

J = r x (— iV — sA) — s- 



(2.12) 



(2.13) 



(2.14) 



J z = s — id/dtp its z-component and J z tj) = mi/). The operator M is written in the Cartesian coordinates as 



M = -r 2 A + XiX 



d 2 



' 3 



dxidxj 



d 2isr ( d d c 2 
2xi- 1 x- y- is - i— 

oxi r — z \ ay ox s 



2c x r 



(2.15) 



Let us consider the generalized MIC-Kepler system in the parabolic coordinates coordinates £, rj € 
[0, oo), ip G [0, 2ir), defined by the formulae 



i'/sin^, z=i(f-r?), 



'tT}CDB(p, y 

In this coordinates the differential elements of length and volume read 
dl 



l 



while the Laplace operator looks like 



A = 



»7 



— (<—)+ — ( — 



1 9 2 

£?7 dip 2 



The substitution 



a i(m— s)ip 

V2n 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



3 



separates the variables in the Schrodinger equation and we arrive at the following system of equations 



v- 



K m? 1_ 1 



£7 mk \ n \ 



477 



$ 2 



0. 
0. 



where /3 - is the separation constant. 

These equations are analogous with the equations of the hydrogen atom in the parabolic coordinates 
Thus, we get 

e i(m-s)(p 



(2.20) 
(2.21) 



where 



1 



r(m, + 1) y 

Here n\ and n 2 are nonnegative integers 

\m-s\+Si + l f3+l 



T(rii + rrii + 1) _<e 21 , . nil 

- e 2 ( ex ) 2 _F(— n.;; m 4 + 1; ex). 



ui = -- 



n 2 = — - 



|m+s| + <y 2 + l /3-1 



(2.22) 



(2.23) 



(2.24) 



2 2e ' 2 2e 

From the last relations, taking into account H2.11|) . we get that the parabolic quantum numbers ri\ and n 2 
are connected with the principal quantum number n as follows 



n = ni + n 2 + 



\m — s\ + \m + s\ 



1. 



(2.25) 



Excluding the energy E from Eqs. (|5.4|l and (|5.5|l . we obtain the additional integral of motion 



X 



a / a 



a ( d 



dr\ \ drj ) dt; \ dt; 



£ — r\ d 2 . £ 2 + rj 2 d 

— IS 

2£?y dip 2 £,rj(£ + rj) dip 



(2.26) 



2c 2 ^ (-1 



with the eigenvalues 



\m — s\ — \m + s\ + 5\ — 5 2 



/? — e m - n 2 + 



and eigenfunctions Vnin 2 m(^, ?7, V! ^2), i-e. 

In Cartesian coordinates, the operator X can be rewritten as 



(2.27) 



(2.28) 



X = z 



d 2 

dx 2 dy 2 



O 2 



d 2 



d__ d_ 

dxdz a dydz r(r — z) \ dy dx 



d_ 
"d~z 



Cl 



C2- 



r(r — z) r[r + z) r(r — 2) r 



(2.29) 



so that it immediately follows that X is connected to the z-componcnt I z of the analog of the Runge-Lenz 
vector 



2 



(-iV — sA) x J - J x (-iV - sA) 



via 



X = 7 2 + ci 



(-2- 



r(r + z) r(r — z) 



(2.30) 
(2.31) 



and coincides with I z when ci — c 2 = 0. 
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3 Bi-Orthogonality of the Radial Wavefunctions 

We shall prove that along with the condition (|2.9[1 the radial wave functions (r; 5i,5 2 ) satisfy the following 
additional orthogonality condition 



Ijj> = I R ( d(r;S 1 ,S 2 )R { n s Ur;5 1 ,S 2 )dr 



5 33' 



( n+ k±h.y2j + 5 1 +5 2 + l 

This new relation is used in the next Section to derive interbasis expansions. It can be proved as follows 

In the integral appearing in (|3.1I) . we substitute explicit expressions <|2.5|) for R^J (r; 5\, S 2 ) and R^J, (r; 
Then, we take the confluent hypergeometric function in (|2.5[) as an finite sum 



(3.1) 



n—j—1 { . . . 1 

( — 71 + 7 + l\p 

pl{2j + Si + S 2 + l)p 
and perform the integration term by term with the help of the formula |22| 



00 

J e- A 7; kx) dx = F( A ^ 1 1) 2 F l 



k 



(3.2) 



(3.3) 



Applying the formula 



2-F1 {a,b;c; 1) 



for the hypergeometric function, we obtain 
T(j+f + S 1 +S 2 + 1) 



l 33 



T(2j + S!+S 2 + 2) 



n—j—l 

£ 

p=0 



T(c)T(c-a-b) 
r(c- a)T(c-b) 



T(n + j + S 1 + S 2 + l)) 



(n-j- l)!(n-i' - l)lT(n + j' + 5i + 5 2 + 1) 



1/2 



(-n + j + l) p (j + f + Sx + 5 2 + l) p r(n - j - p) 



p\(2j + Si + S 2 + 2) p 



r(i'-i- P +i) 



By introducing the formula 



r(z) r(- g + n + i) 

r(«-n) 1 J r(-ar + l) 



(3.4) 



(3.5) 



(3.6) 



into (|3.ti|) . the sum over p can be expressed in terms of the 2 Fi Gauss hypergeometric function of argument 
1. We thus obtain 



1 



j +j' + 5 1 +5 2 + l 
2 



(n-j-l)\T(n+j + S 1 + S 2 + l)) 



(n - j' - l)!r(n + j' + 5 l + 5 2 + l) 
1 



1/2 



(3.7) 



( n + ii±^) 3 r(j-j' + i)r(j'-j + i)' 



Equation (j!H)l then easily follows from (|3~%|) since [T(j - j' + l)T(j' - j + 1)] 1 = 

The result provided by Eq. (|3.1|) generalizes the one for the hydrogen atom [23]. Indeed, orthogonality 
properties similar to (|3.1() hold for the Kepler-Coulomb system and harmonic oscillator in /-dimensional 
spaces (/ > 2) 24 . 
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4 Interbasis Expansion Between Parabolic and Spherical Bases 

The connection between spherical (r, 9, ip) and parabolic (£, rj, tp) coordinates is 

£ = r(l + cos 77 = r(l + cos 0), tp(parabolic) = tp(spherical) . (4-1) 

Now, we can write, for fixed value energy £7„ , the parabolic bound states (|2.22[1 as a coherent quantum 
mixture of the spherical bound states (|2.3|) 



n-1 



^iU(?> »7, ¥>! <*i, £2) = W nm 2 ms (Si, S2) ^ijm (r, 9, ip; 5i,S 2 ) 



(4.2) 



3=m + 



By virtue of Eq. (|4.ip . the left-hand side of 14. 2|) can be rewritten in spherical coordinates. Then, by 
substituting = in the so-obtained equation and by taking into account that 



p(-,« (1)= (EL±l)n 



(4.3) 



we get an equation that depends only on the variable r. Thus, we can use the orthogonality relation (| 3 - 1 ft 
on the quantum number j. This yields 



W^ n2ms (Sx,S 2 ) 



y/(2j + Si + S 2 + l)(j-m+y. EJms Rnni 



where 



and 



Tsnm 
jms 



r(mi + l)r(2j + ft + <5 2 + 2) ~ ni ™ 2 Jms ' 
r (j -m-+Si + l) T(m + m x + i)r(n 2 + m 2 + i)r (71 + j + 5i + <5 2 + 1) 



(n x )!(n 2 )!(n - j - 1)T (j + m_ + <5 2 + 1) T (j + m + + ft + 8 2 + 1) 



-V +mi+5l+52 f 1 (-ni;mi + l;x) F (-n + j + 1; 2j + tfi + <5 2 + 2; a;) dx. 



(4.4) 



(4.5) 



(4.6) 



To calculate the integral K"™*, it is sufficient to write the confluent hypergeometric function F (— m;mi + 1; x) 
as a series, integrate according to (|3.3|) and use the formula l|3.4|l for the summation of the hypergeometric 
function %Fi. We thus obtain 

nni = (n - m+ - l)!r(2j + ft + <5 2 + 2)r (j + m+ + S 1 + S 2 + 1) 
3ms {j-m+)\T(n + 3 +81+82 + 1) 

x 3 ^ 2 

The introduction of l|4.5|l and (|4.7H into l|5.9|l gives 
W j (Si 5?) : 

' r niTi2ms V^l;^ 2 / 



— m, — j + m + ,j + m + + <5i + 8 2 + 1 
mi + 1, — n + m + + 1 



/ (2j + S 1 +S 2 + 1) r(m + mi + l)r(n 2 + m 2 + 1) 
/ (m)!(n 2 )!(n - j - l)!(j - m+)!r (j + m_ + <5 2 + 1) 

{n-m+-l)\ It (j - m_ + ft + 1) T (j + m+ + ft + ft. + 1) 



(4.7) 



r(mi + f) 



x 3 ^ 2 



T(n + j + S 1 +8 2 + l) 

— ni, —j + m+, j + m + + Si + 8 2 + 1 
771 1 + f , — n + 771 _|_ + 1 



(4.8) 
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The next step is to show that the interbasis coefficients (|5.1U|I are indeed a continuation on the real line 
of the Clebsch-Gordan coefficients for the group SU(2). It is known that the Clebsch-Gordan coefficient 



Ca a-b a can written as 



°act;b/3 



(2c+l)(a + a)!(c + 7)! 



(a - a)!(c - 7)! (a + 6 + c + l)!(a + b - c)!(a - 6 + c)!(6 - a + c)! 



1/2 



(4.9) 



X ( _1 J d 7,Q+/3 7= == ; =^= 3-^2 



By using the formula 



3i?2 ^ t,l-N-t 



s,s'-N 



-a — b — c — 1, —a + a, — c + 7 
-a — b + 7, —6 — c + a 



(t + s) N 
(*) 



N 



3F2 



s,t'-s',-N 
t',t + s 



(4.10) 



equation l|4.1U[) can be rewritten in the form 

(2c + 1)(6 - a + c)!(a + <*)!(& + /3)!(c + 7)! 



^,£7 _ 
°aa;b/3 ~~ 



(6 - /3)!(c - 7)! (a + 6 - c)!(a - 6 + c)!(a + 6 + c + 1)! 



1/2 



(4.11) 



(_l)a-^ (q + b-7)! f -a + a,c + 7 + l,- C + 7 
%Q+/3 7(^)!(&-a + 7)! 3 2 I 7-a-6,6- a + 7 +l 

Finally, comparing l|4.12fl and (|5.1U[) . we obtain the representation 

S l + 5 2 rn 1 +m 2 



+ m_+S 2 -l m2 + n 2 -r. 



m.,4- 



1 +n 1 -n 2 



(4.12) 



Equation (|4.12l) proves that the coefficients for the expansion of the parabolic basis in terms of the spherical 
basis are nothing but the analytical continuation, for real values of their arguments, of the SU(2) Clebsch- 
Gordan coefficients. 

The inverse of Eq. (|4.2(l , namely 



n—rriA. — 1 



(4.13) 



is an immediate consequence of the orthonormality property of the SU (2) Clebsch-Gordan coefficients. The 
expansion coefficients in (|4.13[) are thus given by 



1+^2 m 1 +m 2 



, m+\m— s\- 



(4.14) 



and may be expressed in terms of the 3F2 function through i|4.10[l or l|4.12[l . 



5 Prolate Spheroidal Basis 

We now pass to the prolate spheroidal coordinates 

x = "^"VV 2 - !)(! ~ ^ 2 )cos<yj, y = -^VO 2 - !)(! - ^ 2 )sin^, z = — (^+1), (5.1) 

where [i g [0; 00), v G [—1; 1], ip S [0, 27r), and i? g [0; 00). The parameter i? is the interfocus distance, and in 
the limits where R — > and R — ► 00, the prolate spheroidal coordinates give back the spherical coordinates 
and the parabolic coordinates, respectively |21l I27j . 
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The Laplace operator in these coordinates has the form 
4 



i? 2 (/J 2 -1/2) 

After the substitution 



9 2 



i? 2 (/z 2 - 1) (1 - i/ 2 ) dy? 2 



27T 



the variables in the Schrodinger equation l|l.lf> are separated 



d/i 



[f 1 



rfjtx 2(/i+l) 2(/i-l) 



A fl _ 2 ) A _ 1 m 2 

^ ^ > dv 2(1 + j/) 2(1 -i/) 



Rn- 



Rv 



ER 2 
' ~~ 2~ 

ER 2 



(1 - , 2 ) 



= A(i?)Vi 



V>2 = -A(i?)V>2 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



where A(i?) is a separation constant in prolate spheroidal coordinates. By eliminating the the energy E from 
Eas. l|5.4|l and i|5.5|l . we produce the operator 



A 



[i 2 — v 2 



i 2- M 2 -, 2 d 2 (^ + ,) 2 -( M + i)(i + ,) / a 

^ (H 2 - l) (1 - v 2 ) dcp 2 + " + i/) ( M - 1) (1 - i/) v dip;' 
_ 2ci (m + ^) 2 + (m-i)(i-^) + 2c2 (^ + ^) 2 -(m + i)(i + ^) + fi + 1 



(5.6) 



(/i + i/) + i/) 



+ i/) (^i — 1) (1 — /i + ^ ' 



the eigenvalues of which are A(i?) and the eigenfunctions of which are t/Ka'j ^ V 5 )- The significance of the self- 
adjoint operator A can be found by switching to Cartesian coordinates. Passing to the Cartesian coordinates 
in 15.7|l and taking l|2.15|l and (|2.29|) into account, we obtain 



Therefore, 



A = M + RX. 



A ^n<L(M, v,ip\R,5i,5 2 ) = \ q (R)ip$ m (fi, v,(p;R,8 1 ,S 2 ), 



(5.7) 



(5.8) 



where index q labels the eigenvalues of the operator A and varies in the range < q < n — m + — 1. 
We are now ready to deal with the interbasis expansions 

n-l 

v, tp;R,6 1 ,5 2 ) = ]T W nqms (R; S 1 ,S 2 ) V>gn (r, 9, i P ;S 1 ,5 2 ). 

j=m + 



(5.9) 



n— 77i _|_ — 1 
ni=0 



(5.10) 



for the prolate spheroidal basis in terms of the spherical and parabolic bases. (Eq. (|5.9|) was first considered 
by Coulson and Joseph |2E] in the particular case s = 5i = 5 2 = 0.) 

First, we consider Eq. (|5.9|l . Let the operator A act on both sides of l|5.9|) . Then, by using Eqs. (|5.7() . 
(|5.8|) . and l|2.12|l as well as the orthonormality property of the spherical basis, we find that 



, p , , . , S 1 +S 2 \ ( , 6i + 5 2 
\{R) - I J + — o — )\J + — o — + 1 



Ul qms (R;S 1 ,S 2 )=R £ U?; qms (R;Si,S 2 )(X) jr , (5.11) 



8 



where 



Wij> = I (T, 9, <p; S u 6 2 ) Xi>$ m (r, 0, <fi dV. (5.12) 



The calculation of the matrix element (X)jj' can be done by expanding the basis in l|5.12(l in terms of the 
parabolic wavefunctions [see Eq. (|4.13|) ]and by making use of the eigenvalue equation for X [see Eq. <|2.28[l ]. 
This leads to 



n-m+-l , r. i s 



2n + 5 1 + S 2 



ni=0 



njm nj / m' 



Then, by using Eq. (|4.14|l together with the recursion relation [2*5] 



4c 2 (2c+ l)(2c- 1) 



(c + 7)(c - 7) (6 ~ a + c)(a - b + c)(a + b - c + l)(a + b + c + 1) 



1/2 



(5.13) 



(c- 7 - l)(c + 7 - l)(b- a + c- l)(a - b + c- l)(a + b - c + 2) (a + 6 + c) 

4(c- l) 2 (2c- 3)(2c- 1) 



1/2 



x (5.14) 



(a - /3)c(c - 1) - 7 a(a + 1) + 7 b(b + 1) 1>7 1 
2c(c-l) "^^ ' 



and the orthonormality condition 



Er c 7 /V 7 ' _ r r 

«+/3=7 



(5.15) 



we find that (X)jji is given by 
where 



2 ^J'+U 4i A \ | ( TO i +m 2 )(mi -m 2 ) 



2n + 5i + S 2 



{2j + 5 l +5 2 ){2j + 8 l +5 2 + 2) 



A 3 



(j - m+)(J + m + + <5i + £ 2 )(j - m_ + ^i)Q" + m_ + <5 2 )(n - j)(n + j + + 5 2 ) 
(j + (2j + *i + ^2 - l)(2j + & + 5 2 + 1) 



n 1/2 



(5.17) 



Now by introducing (|5.16|) into l|5.11|l . we get the following three-term recursion relation for the coefficient 



U 3 

w nqms 



, , . , $i + S 2 \ ( . , Si + S 2 , A fl(mi + ma) (mi - m 2 ) 

- L? + — o — J + — o — + 11- 



(2j + S 1 +S 2 ){2j + 5 1 +S 2 + 2) 



2R 



2n + S 1 + 5- A 



[Ai£U^ s (R; 5 1 ,5 2 ) + A^U^R; 6 lt S 2 )] = 0. 



(5.18) 



The recursion relation (|5.18|) provides us with a system n — m + linear homogeneous equations which can be 
solved by taking into account the normalization condition 



71—1 



]T \Ui qm8 (R;S 1 ,5 2 )\ =1. 



(5.19) 



]-m + 



Th eigenvalues A g (i?) of the operator A then follow from the vanishing of the determinant for the latter 
system. 
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Second, let us concentrate on the expansion l|5.10|l of the prolate spheroidal basis in terms of the parabolic 
basis. By employing technique similar to the one used for deriving Eq. I|5.11|l . we get 



X q (R) 



2R 



2n + <5i + S 2 



V^ m8 (R; S U S 2 ) = V mms{R\ S 1 ,6 2 )(M) ni< , (5.20) 



n' =0 



where 



(M) nin[ = / ^:> 2m (£, r,, <p; Si,S 2 ) M^ n V n , m (£, n, 9; 81,82) 



dV. 



(5.21) 



The matrix elements (M)„ 1 „/ i can be calculated in the same way as (X)jji except that now we must use the 
relation [25] 



>(c + 1) - a(a + 1) - b(b + 1) - 2a(3] C^. 6>/3 = 



y/(a + a)(a -a + l)(b - (3)(b + + 1) C^_ 1; w + 



(5.22) 



+ y/(a - a )(a + a + l)(b + f3)(b -0+1) C^.^, 
and the orthonormality condition 

a+b 

E C al;bp C lZ';bl3' = <W<W 
c=|"y| 

permit deriving the formula for the matrix element [M) nin i 

(M) nin; = (m + l)(n 2 + ro_) + (n - m + S 2 )( ni + \m-s\ + S 2 ) + ^(81 - S 2 ){5 1 -5 2 -2) + 



(5.23) 



+to_(to + + <5 2 ) 



5 n ' iTll - \/n 2 (ni + l)(ni + \m - s\ + Si + l)(n 2 + \m - s\ + S 2 )S n ^ ni 



+1 



(5.24) 



-y/ni(n 2 + l)(™i + \m- s\ + 5 X + l)(n 2 + \m - s\ + S 2 + l)8 n ^ ni ^. 
Finally, the introduction of ((5.2511 into (|5.2UI) leads to the three-term recursion relation 



(m + l)(n 2 + to-) + (n - ni + <5 2 )(m + \m -s\ + S 2 ) + -(Si - 5 2 )(5i -S 2 -2) + 



+to_(to + + 8 2 ) + 



2R 



2n + 81 + 8 2 



m 1 -m 2 
m - n 2 H ) - Xq(R) 



Vnqms{ R l8l,8 2 ) 



y/n 2 ( ni + l)(m + \m- s\ + Si + l){n 2 + \m - s\ + 5 2 ) V^+l (R;5i,5 2 ) + 
+ V n i( n 2 + l)(ni + \m - s\ + Si + l)(ra 2 + |to - s\ + S 2 + 1) V^qml s i , 82) 



(5.25) 



for the expansion coefficients V^ 1 ms (i?; i5i, <5 2 ). This relation can be iterated by taking account of the nor- 
malization condition 



n— m-f. — 1 

E Kj ms (R;8i,8 2 )\ 2 = i. 

711=0 



(5.26) 



Here again, the eigenvalues X q {R) may be obtained by solving a system of n — m+ linear homogeneous 
equations. 

It should be mentioned, that the formulae ((4.12(1 and (|4. 14|1 and three-term recursion relations 1)5.18(1 and 
((5.18(1 generalize the analogical results for the following systems: 
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• Hydrogen atom [HI EH 131 123 HO EH , when s = 8i = 8 2 = 0. 

• Generalized Kepler-Coulomb system ^7], when s = 0, 8i ^ 62 7^ 0. 

• Hartmann system |34| . when s — 0, 8i = 5 2 =/= 0. 

• Charge-dyon system [S], when s ^ 0, #i = S 2 = 0. 
Finally, it should be noted that the following four limits 

l i ^Uj nqms (R-,S 1 ,8 2 ) = S jq , lim Ul (R;S 1 ,S 2 ) = W^ in ms (S 1 ,S 2 ), (5.27) 

lim V^„.(J2; fc, fc) = S niq , lim V™,,(ii; fc, <5 2 ) = W£ ms (<$i, <5 2 ) (5.28) 
furnish a useful means for checking the calculations presented in Sections 4 and 5. 

Acknowledgements. I would like to thank Dr. Armen Nersessian for useful discussions. The work is 
carried out with the support of ANSEF No: PS81 grant. 

References 

[1] L.G. Mardoyan. The generalized MIC-Kepler system, quant-ph/ 0306168; J. Math. Phys., 44, No. 12 
(2003) (in press). 

[2] D. Zwanziger. Phys. Rev., 176, 1480 (1968). 

[3] H. Mcintosh and A. Cisneros. J. Math. Phys., 11, 896 (1970). 

[4] T. Iwai and Y. Uwano. J. Phys. A. 21, 4083 (1988). A. Nersessian and V. Ter-Antonyan. Mod. Phys. 
Lett. A. 9, 2431 (1994); Mod. Phys. Lett. A., 10, 2633 (1995). 

[5] A. Nersessian, V. Ter-Antonyan, MM. Tsulaia. Mod. Phys. Lett. A. 11, 1605 (1996); A. P. Nersessian, 
V.M. Ter-Antonyan. Physics of Atomic Nuclei, 61, 1756 (1998). 

[6] T.Iwai. J. Gcom. Phys., 7, 507 (1990). L.G. Mardoyan, A.N. Sissakian and V.M. Ter-Antonyan. Physics 
of Atomic Nuclei, 61, 1746 (1998); 

[7] C.N. Yang. J. Math. Phys., 19, 320 (1978). 

[8] L.G. Mardoyan, A.N. Sissakian and V.M. Ter-Antonyan. Int. J. Mod. Phys. A. 12, 237 (1997). 

[9] L.G. Mardoyan, A.N. Sissakian and V.M. Ter-Antonyan. Theor. and Math. Phys., 123, 451 (2000); 
L.G. Mardoyan. Physics of Atomic Nuclei, 65, 1063 (2002). 

[10] V.V. Gritsev, Yu.A. Kurochkin and V.S. Otchik, J. Phys. A. 33, 4903 (2000). 

[11] A. Nersessian and G. Pogosyan. Phys. Rev. A. 63, 020103(R) (2001). 

[12] I. Mladenov and V. Tsanov. J. Phys. A 20, 5865 (1987). 

[13] T. Iwai and Y. Uwano. J. Math. Phys. 27, 1523 (1986). 

[14] A. Inomata, G. Junker, and R. Wilson. Found. Phys. 23, 1073 (1993). 

[15] S. Bellucci and A. Nersessian. Phys. Rev. D, 67, 065013 (2003). 

[16] L.G. Mardoyan, L.S. Petrosyan, and H.A. Sarkisyan. Phys. Rev. A, 68, 014103 (2003). 

[17] M. Kibler, L.G. Mardoyan and G.S. Pogosyan. Int. J. Quan. Chem., 52, 1301 (1994). 



11 



[18] H. Hartmann. Theor. Chim. Acta 24, 201 (1972); H. Hartmann, R. Schuch, and J. Radke. Theor. Chim. 
Acta 42, 1 (1976); H. Hartmann and R. Schuch. Int. J. Quant. Chem. 18, 125(1980). 

[19] J. Fris, V. Mandrosov, Ya.A. Smorodinsky, M. Uhlir, and P. Wintcrnitz. Phys. Lett. 16, 354, 
(1965); P. Winternitz, Ya.A. Smorodinsky, M. Uhlir, and J. Fris. Sov. J. Nucl. Phys. 4, 444, (1967); 
A. A. Makarov, J. A. Smorodinsky, Kh. Valiev, and P. Winternitz. Nuovo Cimento A 52, 1061, (1967). 

[20] N.W. Evans. Phys. Lett. A 147, 483, (1990); Phys. Rev. a 41, 5666, (1990); J. Math. Phys. 31, 600, 
(1990). 

[21] I.V. Komarov, L.I. Ponomarev, and S.Yu. Slavyanov. Spheroidal and Coulomb Spheroidal Functions 
(Nauka, Moscow, 1976). 

[22] L.D. Landau, E.M. Lifshitz. Quantum Mechanics (Pergamon Press, Oxford, 1977). 

[23] A. Erdelyi, W. Magnus, F. Obcrhcttinger, and F. Tricomi. Higher Transcendental Functions (McGraw- 
Hill, New- York, 1953) Vol. I. 

[24] L.G. Mardoyan, G.S. Pogosyan, and V.M. Tcr-Antonyan. Sov. Contemp. Phys., 19, 3 (1984). 

[25] D.A. Varshalovich, A.N. Moskalev, and V.K. Khersonskii. Quantum Theory of Angular Momentum. 
World Scientific, Singapore, (1988). 

[26] W.N. Bailey. Generalized Hypergeometric Series. Cambridge Tracts No. 32, Cambridge, (1935). 

[27] L.G. Mardoyan, G.S. Pogosyan, A.N. Sissakian, and V.M. Ter-Antonyan. J. Phys. A16, 711 (1983). 

[28] C.A. Coulson and A. Joseph. Proc. Phys. Soc, 90, 887 (1967). 

[29] M. Kiblcr, G. Grenet. J. Math. Phys., 21, 422 (1980). 

[30] D. Park. Z. Phys., 159, 155 (1960). 

[31] C.B. Tarter. J. Math. Phys., 11, 3192 (1970). 

[32] M.G. Arutyunyan, G.S Pogosyan, and V.M. Tcr-Antonyan. Sov. Contemp. Phys., 13, 235 (1978). 

[33] L.G. Mardoyan, G.S. Pogosyan, A.N. Sissakian, and V.M. Tcr-Antonyan. Theor. Math. Phys., 64, 171 
(1985). 

[34] I.V. Lutsenko, G.S. Pogosyan, A.N. Sissakian, and V.M. Ter-Antonyan. Theor. Math. Phys., 83, 419 
(1990). 



12 



